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Abstract 

Blazars are a class of active galactic nuclei (AGNs) known for their very rapid variabilty in 
the high energy regions of the electromagnetic spectrum. Despite this known fast variability, 
X-ray observations have generally not revealed variability in blazars with rate doubling or 
halving timescales less than approximately 15 min. Since its launch, the Swift X-ray Tele¬ 
scope has obtained 0.2-10 keV X-ray data on 143 AGNs, including blazars, through intense 
target of opportunity observations that can be analyzed in a multiwavelength context and 
used to model jet parameters, particularly during flare states. We have analyzed this broad 
Swift data set in a search for short timescale variability in blazars that could limit the size 
of the emission region in the blazar jet. While we do find several low-significance possible 
flares with potential indications of rapid variability, we find no strong evidence for rapid (<15 
minutes) doubling or halving times in flares in the soft X-ray energy band for the AGNs analyzed. 


Subject headings: acceleration of particles - galaxies: active - galaxies: jets - gamma rays: general - 
X-rays: general 


1. Introduction 

Analysis of the flux variability of active galactic 
nuclei (AGNs) allows constraints to be put on the 
associated emission mechanisms and jet parame¬ 
ters of the AGN jet. Flux variability in AGNs 
is seen across the entire electromagnetic spectrum 
with the fastest variability seen in the higher en¬ 
ergy regions, in particular the keV to TeV regions. 
Of all AGNs, the strongest and fastest variability 
is often seen in a class of AGNs known as blazars, 
which are AGNs that have a jet pointed in a direc¬ 
tion that is very close to our line of sight. Analysis 
of the flux variability of blazars in the high-energy 
region of the electromagnetic spectrum is a useful 
method for limiting the size of the emission region 
in an AGN, thereby limiting the emission mecha¬ 
nisms. 

Due to light propagation time and geometry 
effects, fast variability can limit the size of the 
emission region R of a blazar by the inequality 


R < (cTS)/{ 1 + z) where T is the variability 
timescale of the flux, S is the Doppler factor, and 
z is cosmological redshift. The Doppler factor 8 
is defined as (5 _1 = y(l — ficosO), where 7 is the 
Lorentz factor in the jet, j3 is the ratio of the jet 
velocity to the speed of light in vacuum, and 0 is 
the angle of the jet with respect to the observer. 
The variability timescale, T, is typically described 
by the flux doubling time, which will be referred 
to as T double or the flux halving time, which will 
be referred to as T^. 

There has been repeated evidence of fast flux 
variability in the TeV 7 -ray region with doubling 
or halving times as short as 224±60 s from PKS 
2155-304 (Aharonian et ah, 2007), 15 minutes 
from Mrk 421 (Gaidos et ah, 1996) and ~ 9 min 
from BL Lac (Arlen et ah, 2012). Results from 
simultaneous TeV and X-ray studies probing vari¬ 
ability on the order of hours have shown evidence 
that X-rays and TeV 7 -rays may originate from 
the same region (e.g. Maraschi et al. 1999). This 
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correlation at larger timescales, in addition to ob¬ 
served fast variability of TeV 7 -rays, has moti¬ 
vated the search for fast variability in blazars at 
X-ray wavelengths. 

In general, X-ray observations have not re¬ 
vealed variability with doubling timescales less 
than approximately 15 min. There has been iso¬ 
lated evidence for much faster variability in X-rays 
from a single blazar, H0323+022 (Feigelson et al., 
1986), with a halving time less than 30 seconds 
and a quoted statistical significance of P < 10~ 13 . 
However, this fast variability in X-rays was a sin¬ 
gle, isolated event that has not been seen again. 
The authors did consider systematic and/or in¬ 
strumental effects, but none could be identified. 

With the launch of the Swift X-Ray Telescope 
(Burrows et al., 2005), X-ray data in the 0.2-10 
keV region has been obtained on many AGNs on 
timescales ranging from seconds to over 8 years. 

We report on our search for fast variability in 
blazar flares. We have searched our large Swift 
database for blazar flares with doubling or halv¬ 
ing times less than 15 minutes. The data taken 
from the Swift database for analysis are outlined 
in Section [2] The analysis techniques, including 
the necessary constraints set on the data in the 
search for fast blazar variability, is outlined in Sec¬ 
tion [3j We provide the results of the analysis in 
Section [4] with discussion of the results in Section 

El 

2. Data 

These data were obtained with the X-Ray Tele¬ 
scope (XRT) on the Swift observatory (Burrows 
et al. 2005), and the reduced light curves are all 
available on a public web site and database de¬ 
scribed by Stroh & Falcone (2013). The Swift- 
XRT is one of three instruments onboard the Swift 
Gamma Ray Burst Explorer (Gehrels et al., 2004). 
In particular, the Sm/f-XRT is an X-ray imag¬ 
ing spectrometer with the ability to obtain data 
for lightcurves with high timing resolution. The 
Swift-XRT data in this study are from two differ¬ 
ent modes, window timing mode (WT) and pho¬ 
ton counting mode (PC). All data are screened to 
ensure that no individual observations, and there¬ 
fore no potential flares, include a switch from one 
mode to the other. See Burrows et al., 2005 for an 
in depth discussion of the XRT modes. 


The data are comprised of 7.43x10® seconds 
of AGN observations from 12544 continuous ob¬ 
servations of the 143 AGN observed at the time 
of analysis. These data were obtained between 
2004 December 17 and 2014 May 28. The data 
were processed using the most up to date version 
of Swift tools at the time of analysis: Swift Soft¬ 
ware version 4.1 and FTOOLS version 6.14 (Black¬ 
burn, 1995) with observations processed using xrt- 
pipeline version 0.12.8. For WT mode data, we 
additionally filtered out the first 150 seconds of 
data after a Swift telescope slew since the tele¬ 
scope pointing was not always settled well enough 
to allow accurate standard WT mode data anal¬ 
ysis. Greater detail on how these data were ob¬ 
tained and processed is available in Stroh & Fal¬ 
cone (2013). 

The number of observations for each object 
range from as many as 1192 observations in the 
case of Mrk 421 to as little as 1 observation in the 
case of PKS 0244-470. The average continuous 
observation time is about 592 s, with a maximum 
observation time of 2565 s and a minimum obser¬ 
vation time of 21 s. 

3. Analysis 

The goal of the analysis is to determine if fast 
blazar variability with a rate doubling or halv¬ 
ing time less than 15 minutes exists within these 
data. The data were analyzed on an observation 
by observation basis to test for this fast variability. 
Various constraints were needed to: (1) define the 
data set of observations that would have enough 
counts and a long enough duration to provide the 
potential to find significant flaring, and ( 2 ) define 
the observations within the aforementioned data 
set that would have variability, and in particular, 
the potential for flux doubling or halving. The 
necessary constraints to define our data set in¬ 
clude: a maximum count rate of at least 0.1 cts/s, 
at least three data points in each observation, and 
a length of observation of at least 5 minutes. The 
constraints that define a possible flare include: re¬ 
stricting the maximum count rate to be at least 1.8 
times greater than the minimum count rate during 
the observation, 3 consecutive data points either 
increasing or decreasing in count rate, and a re¬ 
duced x 2 value from a straight line fit greater than 
3.0. In addition to these pre-defined constraints, 
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Count Rate (0.3-10.0 keV) (cts s' 1 ) 


0954+658, Obs ID: 00035381001 1ES 2344+514, Obs ID: 00035031050 3C 279, Obs ID: 00035019131 



53920.440 53920.445 53920.450 53920.455 55199.415 55199.420 55199.425 55199.430 56657.475 56657.480 56657.485 56657.490 

Time (MJD) Time (MJD) Time (MJD) 


3C 279, Obs ID: 00053516001 3C 345, Obs ID: 00036390005 3C 454.3, Obs ID: 00031493027 



53926.200 53926.205 53926.210 53926.215 55066.805 55066.810 55066.815 55066.820 56567.775 56567.780 56567.785 56567.790 

Time (MJD) Time (MJD) Time (MJD) 



PKS 1424+240, Obs ID: 00041539012 PKS 1510-089, Obs ID: 00031173069 RX J0324.6+3410, Obs ID: 00036533033 



56395.310 56395.315 56395.320 56395.325 55679.855 55679.860 55679.865 55679.870 56307.855 56307.860 56307.865 56307.870 

Time (MJD) Time (MJD) Time (MJD) 


S5 0716+714, Obs ID: 00035009047 S5 0716+714, Obs ID: 00035009048 W Com, Obs ID: 00035018040 



55179.700 55179.705 55179.710 55179.715 55180.295 55180.300 55180.305 55180.310 55356.735 55356.740 55356.745 55356.750 

Time (MJD) Time (MJD) Time (MJD) 


Fig. 1. The light curves of the 15 observations flagged for a potential quick time X-ray flare are shown. 
The observation ID specified in the title of each plot is the observation ID assigned by Swift. 
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manual inspection was required for the observa¬ 
tions that passed the constraints to ensure that 
flux variations from systematic effects (e.g. mode 
switching) were not flagged as possible flares. The 
necessity of each of these constraints is discussed 
in more detail below. 

3.1. Constraints That Define Data Set 

3.1.1. Max Count Rate 

A maximum count rate of at least 0.1 cts/s dur¬ 
ing an observation is necessary to allow for a sta¬ 
tistically significant flare to occur in the timescales 
being probed. Taking into account light curve 
error bars at this count rate and applying y 2 
statistics, a statistically significant flare cannot be 
detected at rates lower than 0.1 cts/s over the 
timescales of concern. Therefore, if the maximum 
count rate of an observation is not at least 0.1 
cts/s then it would be impossible to detect a fast 
flare from that observation, thus we eliminate such 
observations from the data sample. 

3.1.2. Data Point Requirement 

The requirement of at least three data points 
within an observation as a means of defining our 
data set is closely related to the discussion in Sec¬ 
tion 13.2.21 We cannot include observations that 
have less than three data points as part of our 
data set because they do not have the ability to 
have three increasing or decreasing data points 
and therefore would not possibly be able to be 
flagged as a possible flare. 

3.1.3. Orbit Length 

We also require observations to be at least 5 
minutes in length. This is to allow enough time 
for a significant flare to occur within a typical ob¬ 
servation. Our data set also has max observation 
time of 2565 seconds, which implicitly sets an up¬ 
per limit to the length of the doubling or halving 
time observed of a potential flare, though no ex¬ 
plicit cut was set. 

3.2. Constraints That Define Flaring 

3.2.1. Doubling/Halving Factor 

We require the maximum count rate to be 1.8 
times greater than the minimum count rate to en¬ 
sure that any potential flares at least approach a 


flux doubling or halving since the purpose of the 
analysis is to search for AGN flares with fast dou¬ 
bling or halving timescales. 

3.2.2. Consecutive Points 

We require at least 3 consecutive data points 
to be either increasing or decreasing in count rate 
to minimize random noise in the count rate and 
guarantee a possible steady increase or decrease 
in the flare. 

3.2.3. Reduced x 2 Value 

Finally, a reduced y 2 value from a straight line 
fit of the weighted mean of the count rate greater 
than 3.0 is necessary to ensure sufficient variabil¬ 
ity in the count rate of the source. This constraint 
ensures that the selected observation is not con¬ 
sistent with a quiescent state of the AGN. This 
reduced y 2 value was also used in determining the 
p-value of the null hypothesis of a straight line fit. 

3.3. Manual Analysis 

Manual inspection was required to ensure that 
all observations flagged as possible flares did not 
have any issues that could not be flagged systemat¬ 
ically. One of these issues includes the Swift-X RT 
mode switching from PC to WT or vice versa, 
which often caused an abrupt break in the light 
curve and spurious points to occur in the data as 
shown in light curve on the top of Figure [2] This 
issue would occasionally cause observations to be 
flagged as a possible flare and would have to be 
discarded manually. Other observations that were 
processed out through manual inspection had a 
single spurious data point responsible for the iden¬ 
tification of the observation as a potential flare, 
with the single point occurring either at the be¬ 
ginning or end of an observation with no flaring 
evidence throughout the rest of the light curve. 
These observations (e.g. Figure 0 bottom) also 
had to be discarded from the data set. 

4. Results 

Of the initial 12544 AGN observations in the 
large data set, 8606 AGNs observations had the 
possibility of having a detectable quick time flare 
based on the count-rate constraints and the a pri¬ 
ori assumptions that define our data set as dis¬ 
cussed in Section EHH Of the remaining 8606 
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Mrk 421, Obs ID: 00030352205 



55242.010 55242.015 55242.020 55242.025 

Time (MJD) 


1156+295, Obs ID:00036381002 



Fig. 2.— The light curve from Mrk 421 (top) 
shows a sharp break in the light curve that results 
in spurious data points as the result of the Swift - 
XRT switching observing modes from PC mode to 
WT mode. The light curve of 1156+295 (bottom) 
shows a single spurious data point that was po¬ 
tentially a result of fluctuations at the beginning 
of an observation that could not be considered a 
potential quick time flare. 

AGN observations, 31 observations passed the 
constraints that defined a possible flare as dis¬ 
cussed in Section 13.21 Once the manual analy¬ 
sis was done, 21 observations remained that did 
not have any mode switching issues or were re¬ 
liant on a single spurious data point, as discussed 
in Section 1X31 Of the remaining 21 observations, 
six were flagged as being clearly due to random 
fluctuations, without any clearly discernible flare, 
and were therefore not considered to be poten¬ 
tial flares, and one observation from PKS 0537-441 
showed signs of two flaring episodes within a single 
observation. Therefore, 16 potential flares from 15 
observations and 12 different AGN remained. De¬ 
tails from these 16 potential flares are shown in 
Tabled and each light curve is shown in the mo¬ 


saic of Figure [I] 

We used two methods to estimate the signif¬ 
icance of the quick time soft X-ray flares; one 
method compares the light curve data to an as¬ 
sumption that the X-ray rate was constant with 
a value equal to the mean rate during the obser¬ 
vation, and a second method compares the light 
curve data to an assumption that the X-ray rate 
was constant with a value equal to a best guess 
at the quiescent source rate near the time of the 
potential flare. In the first method, the chance 
probability p-value was determined based on the 
null hypothesis probability of a straight line fit of 
the weighted mean and using the reduced y 2 value 
that was found as discussed in Section 13.2.31 In 
the second method the chance probability p-value 
was determined based on the null hypothesis prob¬ 
ability of the light curve data being described by 
a constant, straight-line fit to the minimum rate 
surrounding each potential flare, which was used 
as a best-guess estimate of the ’’quiescent” state 
of emission surrounding the potential flare. We 
computed the reduced y 2 value of the data dur¬ 
ing the potential flare relative to this ’’quiescent” 
constant rate light curve and then determined the 
p-value. These values are shown in Table |T] un¬ 
der ’’Weighted Mean p-value” and ’’Quiescent p- 
value” respectively. 

Once the p-value was determined using the two 
different methods, we had to take into account a 
trials factor based on the number of observations 
in the data set that had the potential to provide a 
detection of a rapid soft X-ray flare. The number 
for the trials factor was determined by counting 
the number of observations that passed the set 
constraints as discussed in Section no and Sec¬ 
tion 13.31 as well as the visual inspection described 
above. Therefore, there were 16 trials that were 
taken into account, and the resultant p-values for 
each method are shown in Table HI in the ” Post¬ 
trials” column next to each respective method. 

In Table[l]we also include the reduced y 2 values 
that were found on the section of each observation 
that was deemed part of the potential quick time 
flare as discussed earlier for the weighted mean and 
quiescent straight line fits. We also include the 
doubling times and halving times when applicable 
to each flagged observation. These times were de¬ 
termined by fitting a straight line to the sections 
of each observation deemed to be part of a pos- 
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RX J0324.6+3410, Obs ID: 00036533033 


0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 

56307.855 56307.860 56307.865 56307.870 

Time (MJD) 



1 ES 2344+514, Obs ID: 00035031050 



55199.415 55199.420 55199.425 55199.430 

Time (MJD) 


Fig. 3.- These plots show the selected section of 
the observations from RX J0324.6+3410 on 2013 
Jan 15 (top) and 1ES 2344+514 on 2012 January 
3 (bottom) that were fit with either increasing or 
decreasing lines to estimate the doubling or halv¬ 
ing time. 


sible flare and calculating how much time would 
be required for half of the maximum count rate 
to double or for the maximum count rate to halve 
when appropriate. An example of this fitting can 
be seen in Figure [3] We also include the date of 
observation of each potential flare as well as the 
start time of the observation containing the possi¬ 
ble flare as a reference to compare the data to the 
light curves available on the Swift-XKT monitor¬ 
ing site (Stroh & Falcone, 2013). 

5. Discussion 

Analysis of our large Swift-XKT database in¬ 
dicates that our data set does not show strongly 
significant evidence for soft X-ray flaring on the 
timescales probed, i.e. less than 15 minutes. 
Of the low-significance potential flares that were 
found and outlined in Table |T] and Figure [Q we 


find potential doubling times in 6 different AGNs 
between 2.2 min and 5.7 min as well as halving 
times in 10 different AGNs between 2.2 min and 
9.1 min. All of the potential flares in Tabled] come 
from the blazar class (FSRQs or BL Lacs), which 
are more likely to exhibit rapid flaring due to their 
aligned relativistic jets. However, none of these 
potential flares represents a truly significant de¬ 
tection after accounting for trials, and therefore, 
no strong statements can be made about the exis¬ 
tence, or nonexistence, of short timescale flares in 
the data set. 

The potential flare from the RX J0324.6+3410 
observation is shown in more detail in Figure [3] 
(top). This potential flare has a doubling time, 
calculated as discussed in Section [4] of 5.7 min¬ 
utes. If the variability is truly from flaring, the 
potential emission region size can be calculated 
by applying the inequality R < (cT5)/(l + z), 
which was discussed in Section [T] The doubling 
time variabilitiy as well as z=0.061 (Linford et al. 
2012 ) equate to an upper limit on the emission re¬ 
gion size of 9.7x 10 12 d cm, where S is the unknown 
doppler factor that was discussed in Section [T] 

If the potential flares outlined in Table [1] were 
indeed quick time soft X-ray flares, these would 
be the fastest doubling and halving times for X- 
ray AGN flares ever observed, with the exception 
of the single, isolated event observed by Feigelson 
et al. (1986). As stated previously, variability in 
AGNs at timescales approximately this short had 
previously only been seen repeatedly with statis¬ 
tical significance in the TeV region (Aharonian et 
al., 2007; Gaidos et al., 1996, Albert et al., 2007). 

Simultaneous multiwavelength studies of blazars 
in TeV 7 -rays and X-rays are vital to determine 
if very fast X-ray variability is coincident with 
very fast TeV variaiblity as has been seen previ¬ 
ously on longer timescales (Maraschi et al., 1999). 
Further X-ray observations of these sources with 
more sensitive instruments for longer continuous 
integrations to determine flare timescales will al¬ 
low us to perform better searches for fast X-ray 
variability in the future. This will enable more 
complete studies of the connection between X-ray 
and TeV emission regions in blazar jets. 
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Table 1: The 16 potential quick time soft X-ray flares flagged for more analysis. 


Source 

Date 

Obs. Time 

Weighted 

Mean 

Xred 

Quiescent 

X-red 

T double 

Ti/2 

Weighted 

Mean 

p-value 

Post-Trials 

Quiescent 

p-value 

Post-Trials 

TXS 0954+658 

04-Jul-06 

53920.442 

3.23 

4.38 


9.1 

3.96xl0" 2 

6.34xl0 -1 

1.25xl0” 2 

1.99xl0 _1 

1ES 2344+514 

03-Jan-10 

55199.419 

4.26 

6.70 

- 

2.8 

5.13xl0" 3 

8 .20xl0~ 2 

1.61xl0~ 4 

2.58xl0 -3 

3C 279 

31-Dec-13 

56657.480 

4.45 

11.96 

3.4 

- 

3.95xl0 -3 

6.32xl0~ 2 

7.91xl0 -8 

1.27xl0“ 6 

3C 279 

lO-Jul-06 

53926.208 

3.05 

4.69 

- 

3.9 

2.73xl0~ 2 

4.36X10' 1 

2.81xl0” 3 

4.49xl0" 2 

3C 345 

23-Aug-09 

55066.809 

3.19 

4.29 

- 

7.0 

4.10xl0~ 2 

6.56xl0 _1 

1.37xl0 -2 

2.19xl0 _1 

3C 454.3 

02-Oct-13 

56567.783 

3.63 

11.52 

4.2 

- 

2 .66xl0~ 2 

4.26xl0 _1 

9.90xl0~ 6 

1.58xl0" 4 

PKS 0235+164 

06-Feb-07 

54137.292 

5.95 

7.79 

2.5 

- 

2.60xl0 -3 

4.16xl0” 2 

4.13xl0" 4 

6.61xl0 -3 

PICS 0537-441, 1 

08-0ct-08 

54747.866 

4.37 

6.10 

- 

3.6 

1.27xl0~ 2 

2.03X10" 1 

2.25xl0~ 3 

3.60xl0 -2 

PKS 0537-441, 2 

08-0ct-08 

54747.869 

5.05 

6.76 

3.2 

- 

6.41xl0 -3 

1 .02xl0 _1 

1.16xl0" 3 

1.85xl0 -2 

PKS 0537-441 

26-Jan-05 

53396.605 

5.68 

8.08 

2.2 

- 

3.40xl0~ 3 

5.44xl0 -2 

3.11xl0~ 4 

4.97xl0 -3 

PKS 1424+240 

13-Apr-13 

56395.315 

3.50 

5.45 

- 

4.9 

1.47xl0 -2 

2.35xl0 _1 

9.68xl0 -4 

1.55xl0 -2 

PKS 1510-089 

28-Apr-11 

55679.860 

5.94 

6.96 

- 

5.8 

2.63xl0” 3 

4.21 x 10 -2 

9.46xl0 -4 

1.51xl0 -2 

RX J0324.6+3410 

15-Jan-13 

56307.859 

3.91 

9.46 

5.7 

- 

3.57xl0~ 3 

5.71xl0 -2 

1 .21xl0“ 7 

1.93xl0" 6 

S5 0716+714 

14-Dec-09 

55179.704 

4.05 

6.62 

- 

2.2 

1.74xl0~ 2 

2.78xl0 _1 

1.33xl0” 3 

2.13xl0 -2 

S5 0716+714 

15-Dec-09 

55180.299 

5.51 

10.98 

- 

3.2 

4.05xl0~ 3 

6.48xl0 -2 

1.71xl0” 5 

2.73xl0" 4 

W Com 

09-Jun-10 

55356.738 

3.78 

5.36 

- 

8.1 

2.28xl0~ 2 

3.65xl0 _1 

4.68xl0” 3 

7.49xl0 -2 


Notes.—Col. (1): Object name referred to on Swift-XRT monitoring site. Col. (2): Date of observation. Col. (3): Start time of observation 
in MJD. Col. (4): Reduced % 2 of weighted mean straight line fit. Col. (5): Reduced % 2 of simulated quiescent line fit. Col. (6): Estimated 
doubling time of potential flare in minutes. Col. (7): Estimated halving time of potential flare in minutes. Col. (8): Estimated p-value 
significance of weighted mean straight line fit before trials factor. Col. (9): Estimated p-value significance of weighted mean straight line 
fit after trials factor. Col. (10): Estimated p-value significance of quiescent line fit before trials factor. Col. (11): Estimated p-value 
significance of quiescent line fit after trials factor. 



